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ABSTRACT 

We consider all bursts with known redshift and vF v peak energy, E°^ k . For 
a good fraction of them an estimate of the jet opening angle is available from 
the achromatic break of their afterglow light curve. This allows the derivation 
of the collimation-corrected energy of the bursts, E 1 . The distribution of the 
values of E y is more spread with respect to previous findings, covering about 
two orders of magnitude. We find a surprisingly tight correlation between E 1 
and the source frame E peak : E°^ k (l + z) oc E®- 7 . This correlation can shed light 
on the still uncertain radiation processes for the prompt GRB emission. More 
importantly, if the small scatter of this newly found correlation will be confirmed 
by forthcoming data, it will be possible to use it for cosmological purposes. 

Subject headings: Gamma Rays: bursts — Radiative processes: non-thermal 



The possibility that GRB fireballs are collimated was first proposed for GRB 970508 
(Waxman et al. 1998) and subsequently invoked for GRB 990123 as a possible explanation for 
its extraordinarily large isotropic energy (Fruchter et al. 1999). The observational evidence 
supporting this scenario is the achromatic break of the afterglow light curve which declines 
more steeply than in the spherical case (e.g. Rhoads 1997, Sari, Piran & Halpern 1999). 
Under the simplifying assumption of a constant circum-burst density medium of number 
density n, a fireball emitting a fraction 7/ 7 of its kinetic energy in the prompt 7-ray phase 
would show a break in its afterglow light curve when its bulk Lorentz factor T becomes of 
the order of T ~ 1/8 with 9 given by (Sari et al. 1999) 1 : 
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1 The notation Q = 10 X Q X is adopted, with cgs units. 
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where z is the redshift, tj et: d is the break time in days, and E ljiso is the energy in 7-rays 
calculated assuming that the emission is isotropic. The collimation-corrected energy is 
E 1 = (l- cosd)E~ f , lso . 

Frail et al. (2001 - F01 thereafter) considering a sample of 15 bursts with redshift and 
estimate of 6 (including 5 lower/upper limits) found the remarkable result of a clustering 
of E-y around 5 x 10 50 erg. This was also confirmed independently by Panaitescu & Kumar 
(2001). More recently, Bloom et al. (2003 - B03 hereafter) found that the distribution of E 1 
for a larger sample of 24 bursts (including 8 lower/upper limits) clusters around 1.3 x 10 51 
erg, emphasizing at the same time the presence of several outliers with a sizeably smaller 
energetics. These results suggest that GRBs are characterized by a universal energy reservoir, 
despite the very large range of isotropic equivalent energetics. 

A correlation between E ltiso and the source frame peak energy E peak was discovered 
by Amati et al. (2002) (see also Lloyd-Ronning et al. 2000, Lamb et al. 2003; Sakamoto 
et al. 2004). An analogous correlation holds between the peak luminosity L pea k^ so and 
Ep ea k (Schaefer 2003; Yonetoku et al. 2003). Possible interpretations of these correlations 
(Schaefer 2003; Eichler & Levinson 2004, Liang et al. 2004) are now under intense discussion. 
It seems therefore that the local energy content of the jet plays a crucial role in determining 
the typical photon energy E pea k. In this paper we show that the tightest correlation is instead 
between E pea k and E y . 



2. Sample selection 

We considered all the 40 GRBs with measured redshift up to June 2004 2 . 

The prompt emission spectrum of GRBs can be described by a phenomenological spec- 
tral model (i.e. the Band function, Band et al. 1993) composed by low and high energy 
powerlaws, with photon spectral indices a and (3 respectively, smoothly connected by an 
exponential cutoff with a characteristic energy E break . In the vF v representation this model 
predicts a peak energy E peak = (a + 2)Ei )rea k- For all the bursts with measured redshift 
we searched in the literature any information about their prompt emission spectrum. All 
bursts with firm redshift measurement and published E pea k were included in our final sample 
reported in Tab. 1. This sample contains 29 GRBs detected by different instruments (Col. 1) 
and distributed in redshift (Col. 2) between £=0.0085 (GRB 980425, the second nearest is the 



2 A continuously updated collection of GRBs with links to the relative GCN communications can be found 
at http://www.mpe.mpg.de/jcg/grbgen.html 
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X-ray Flash 020903 with z=0.25, Soderberg et al. 2004) and z=4.5 (GRB 000131, Andersen 
et al. 2000). Col. 3 to Col. 8 report the parameters of the spectrum integrated over the 
burst duration (Col. 6). The average low (Col. 3) and high (Col. 4) energy spectral indices 
of our sample are (a) = —1.05 ± 0.44 and {(3) = —2.28 ± 0.25, which are consistent with the 
typical values found from the spectral analysis of bright long bursts (see, e.g. Preece et al. 
2000). All the GRBs for which a redshift has been measured belong to the population of 
long bursts (e.g. Kouveliotou et al. 1993) and the average duration of the sample reported 
in Tab. 1 is 88 sec. The fluences (although not uniformly integrated over the same energy 
range) are representative of the typical GRB values. For the GRB detected by BeppoSAX 
and studied by Amati et al. (2002) the spectral parameters are derived combining the Wide 
Field Camera (WFC) and the Gamma Ray Burst Monitor (GRBM) instruments on-board 
BeppoSAX, while the fluence is given by the analysis of the GRBM data only. There is 
a small difference between the 7-ray fluences derived fitting the GRBM data only or the 
GRBM+WFC data. Therefore, using the spectral indices and the normalization obtainable 
from the fluences given in Tab. 1 does not lead to the correct value of E lt i S0 reported in Tab. 
2 (for BeppoSAX bursts only). The latter is taken directly from Amati et al. (2002), i.e. 
from the combined fit. Note also that for bright bursts the GRBM instrument is sensitive 
over a larger energy range (extending towards higher energy) than the nominal 40-700 keV 
interval. 

Eight GRBs (Tab. 3) were excluded from this sample either because their redshift was 
uncertain (2 GRBs) or because their peak energy was not found in the literature (6 GRBs). 
They are discussed in Sect. 4.1 and their consistency with the conclusions drawn for the 
main sample is checked under reasonable assumptions for the lacking parameters. 

From the original sample of 40 GRBs with known redshift only 3 cases were completely 
excluded from the present analysis: GRB 020124 (Ricker G. et al. 2002), whose spectral 
properties are reported in Barraud et al. (2003), presents a peak energy of ~ 1 MeV with 
an uncertainty of an order of magnitude; GRB 030323 (Graziani et al. 2003) for which no 
peak energy or fluence was reported in the literature and GRB 031203 (Gotz et al. 2003) 
without a (yet) published prompt phase spectrum. 

Considering all the bursts presented in this work (Tab. 1 and Tab. 3) our sample contains 
37 GRBs. 
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3. Isotropic energy and collimation correction 



The source frame isotropic equivalent energy of a burst can be derived from its frequency 
and time integrated flux. If the jet opening angle is known then we can calculate the 
collimation corrected energy. These quantities for the sample of selected bursts are reported 



Similarly to Bloom et al. (2001, 2003) and Amati et al. (2002) we derived for every 
burst reported in Tab. 1 the source frame "bolometric" isotropic energy E^^o integrating 
the best fit time-integrated model spectrum N(E) [phot cur 2 keV -1 ] over the energy range 
1 keV-10 MeV. The integration over such a large band also required the proper correction 
for the band-redshift effect so that: 



where E is in keV and Di is the source luminosity distance (we adopt f2 m =0.3, Q\=0.7, 
H =70 km s -1 Mpc" 1 ). Note that Amati et al. 2002 assume for the BeppoSAX bursts a 
value of H = 65 km s' 1 Mpc -1 in deriving E 1;iso . For these bursts we corrected their E 1;iso 
(col. 7 in Tab. 2) for our different H . 

For 24/29 GRBs reported in Tab. 1 a jet break time is known and the collimation 
angle 6 can be derived (Eq. 1). For these we used the isotropic energy estimated above 
with a constant energy conversion efficiency ?7 7 =20% (see also F01). Only in few cases the 
circum-burst density was measured (from broad band modeling of the afterglow emission, e.g. 
Panaitescu & Kumar 2001, 2002) and used in deriving 8. For those bursts with unknown n 
we assumed the median value n ~ 3 cm -3 of the distribution of the measured densities which 
extends roughly between 1 and 10 cm -3 (Frail et al. 2000a, Yost et al. 2002, Panaitescu & 
Kumar 2002, Frail et al. 2003, Harrison et al. 2001, Schaefer et al. 2003). This density was 
also assumed for GRB 990123 although Panaitescu & Kumar (2001), fitting simultaneously 
the afterglow light curves in different bands, found for this burst and for GRB 980703 
n ~ 10~ 3 cm -3 . Studying the radio properties of GRB 980703, Frail et al. (2003) found 
n ~ 30 cm -3 , and discussed the possible reasons for the discrepancy with the value quoted 
above (see also Berger et al. 2004a), which could be applied also to GRB 990123 (see also 
Nakar & Piran 2004). Similar to F01 and B03 we have then treated GRB 990123 as a burst 
with unknown circum-bursts density, and therefore used n = 3 cm -3 (note that F01 used 
n = 0.1 cm" 3 , while B03 used n = 10 cm -3 ). 

The jet break time (collected in most cases from B03) is reported in Tab. 2 (Col. 1) 
for all the bursts (16) for which a direct measure of this parameter was possible from the 
broad band modeling of the afterglow light curve. For 7 GRBs only upper/lower limits on 



in Tab. 2. 
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tj e t were estimated from the available data. These limits determine the upper /lower limits 
on the collimation-corrected energy reported in Tab. 1 (Col. 7). On the other hand there are 
a handful of GRBs whose spectrum is best fit by a Band function with high energy spectral 
index f3 >-2. This indicates that E pea k is outside (above) the fitting energy range. These 
cases are reported as lower limits in energy (Tab. 2, Col. 8). Only for XRF 030723 the limit 
on the peak energy and on the isotropic energy is due to the limit on its redshift z <2.1 
(Fynbo et al. 2004). 

The jet break time, the derived jet opening angle, the source frame isotropic energy 
E lt i SO , collimation-corrected energy E 1 and peak energy E pea k are reported in Tab. 2 for the 
sample of bursts in Tab. 1. Tab. 4 reports these quantities for the few GRBs of Tab. 3 with 
uncertain parameters. We note that, due to the larger energy range for the calculation of 
the isotropic energy and to the slightly different circum-burst assumed density, the values of 
9 and E 1 reported in Tab. 2 are different with respect to those reported by F01 and B03 for 
the common GRBs. 

The spectral parameters of several bursts reported in Tab. 1 (and Tab. 2) were pub- 
lished without the relative errors. For these parameters we assumed the average errors (in 
square parenthesis) obtained from those bursts with published parameters' errors (in round 
parenthesis). Furthermore, we do not know the correlations among the parameters entering 
in the calculation of the jet opening angle and of the collimation corrected energy. For this 
reason (see also B03) the simplest approach is to propagate the errors in their calculations 
assuming that these parameters are uncorrelated. In all the tables the errors obtained with 
this method are reported in square parenthesis and the resulting average uncertainty on the 
three foundamental quantities E pea k, 9 and E 7 are 20%, 32% and 20%, respectively. 



In Fig. 1 we report (filled symbols) the source frame peak energy E pea k(l + z) versus the 
collimation corrected energy E 1 for the 24 GRBs with redshift and 9. We find that these 
two parameters are highly correlated with a Spearman's rank correlation coefficient r s =0.88 
and a null hypothesis probability of being uncorrelated of 2.7xl0~ 8 . If we exclude from the 
correlation all the bursts (8) for which either the computed E peak and/or E 7 are upper/lower 
limits we derive a correlation coefficient of r s =0.94 (with a probability of 1.4xl0~ 7 ) and the 
best fit powerlaw model, obtained accounting for errors on both coordinates (routine fitexy 
of Press et al. 1999), is (solid line in Fig. 1) 
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0.706±0.047 



keV 



(3) 
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Note that fitting E y versus E peak with this method gives E 1 oc E^ k 6±om , which is exactly 
equivalent to Eq. 3. Fig. 1 also shows (open symbols), for all the 30 GRBs present in Tab. 1, 
the peak energy and the isotropic equivalent energy E ltiso . This correlation presents a larger 
spread with respect to previous findings (Amati et al. 2002, Lamb et al. 2003) due to the 
larger number of objects included in the sample. The correlation coefficient results r s =0.803 
(with a probability of 7.6xl0~ 7 ) and excluding the upper/lower limits on E pea k (3 GRBs) 
the best fit is E peak = 258(£ 7iiso /1.2 x io 53 er#) a40±a05 [dashed line in Fig. 1). Note that, 
especially for E peak ~ 500-600 keV, there is a relatively large scatter of E^ tiso values around 
the fitting line. For the very same bursts the scatter of E 1 around the fitting E peak — E^ 
line is much reduced. This means that it is possible that bursts with the same and E peak 
have different jet opening angle. 

Of all the bursts included in Tab. 1 only GRB 980425 represents an outlier for the 
E peak -E 1>iso correlation, and for this burst the angle is unknown. Instead XRF030723 and 
XRF020903 are consistent with the extrapolation of the above correlation at very low peak 
energies (see also Lamb et al. 2003). For XRF030723 there are indications of a break in 
the lightcurve between 30 and 50 hours (Dullighan et al. 2003), while for XRF020903 a 
jet break time around 3 days (corresponding to ~ 25 days), would be required to bring it 
on the E peakr E^ correlation. The 8.46 GHz data of XRF020903 presented in Soderberg et 
al. (2004) show a relatively fast time decay (decay index < —1.5), after the source became 
transparent (i.e. ~25 days after the trigger). However, Soderberg et al. (2004) claim that 
there is no sign of jet break at 4.86 GHz, where the flux increases up to ~30 days after the 
trigger. If this is the case, then XRF020903 would be an outlier. 

The correlation that we find between the peak energy and the collimation corrected 
energy is also extremely narrow and we note that except for one burst (GRB 000911) all 
the upper/lower limits are still consistent with this correlation. In Fig. 2(a) we report 
the distribution of the isotropic energy which, similarly to what reported by many authors 
(Bloom et al. 2001, Frail et al. 2001, Berger et al. 2003), is wide spread over almost 3 
orders of magnitude between 6.8xl0 51 erg and 2.8xl0 54 erg. Applying the correction for the 
collimation angle, this distribution clusters around a typical value of ~10 51 erg (Fig. 2(b)). 
The gaussian fit to the logarithmic distribution of E 1 results in a central value (log(i? 7 ))=50.8 
with a standard deviation of 0.6. Introducing the correlation of E 1 with E peak we can derive 
a measure of the dispersion of the GRBs around the correlation E peak oc E®- 7 in the log-log 
plane. The distribution of the dispersion measure is reported in Fig. 2(c) and its average 
value is 0.041 ± 0.015 with a maximum logarithmic dispersion of 0.25. We note that this 
distribution represents the effective tight correlation between E peak and E 1 . 
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4.1. Consistency checks 

From the sample of GRBs with measured redshift we excluded 8 cases with uncertain 
redshift and/or peak energy. These GRBs are reported (with the same meaning of the 
columns) in Tab. 3. In this section we discuss their consistency with the above correlation 
under reasonable assumptions on their lacking parameters. 

GRB 970508 was detected both by BeppoSXX and BATSE. According to the analysis 
of BeppoSAX data (Amati et al 2002) its peak energy should be 145±43 keV whereas the 
analysis of the BATSE spectrum (Jimenez et al. 2001) reports an unconstrained peak energy 
(due to a high powerlaw spectral index (5 >-2). The latter determines only a lower limit 
of E peak >1503 keV (rest frame). A more difficult situation was encountered for 5 GRBs 
(991208, 000210, 000301C, 000418, 000926) with unpublished spectral parameters. Except 
for GRB991208, the other events were detected by one of the satellites of the Interplanetary 
Network. These instruments are most likely sensitive to GRBs with a spectrum similar to 
that of bursts detected by BATSE. For this reasons, having only a fluence, we assumed a 
typical Band function with a and (3 equal to the average values of these parameters for 
the bursts presented in Tab. 1 (see also sec. 2) and with a break energy free to vary in the 
range 100-500 keV which is a conservative estimate of the typical width of the distribution 
of this spectral parameter for BATSE bright bursts (Preece et al. 2000). With the above 
assumptions we iteratively calculated E Jjiso , E peak and E 1 for these bursts and report their 
interval of variation in Fig. 3. In Tab. 4 we report for reference their peak energy, isotropic 
equivalent energy and collimation corrected energy corresponding to the assumption of the 
above average spectrum with break energy equal to the central value of its allowed range of 
variation. 

Finally there are two GRBs (980326 and 980329) with uncertain redshift measurements. 
There are some indications (Bloom et al. 1999) that 980326 lies at z <1.5 (this limiting value 
was assumed by B03) and that 980329 might be between redshift 2 and 3.9. For these two 
events we calculated the region of the E peak -E 1 ^ SO) E 1 plane were they lie varying their z 
within these ranges. They are represented by the connected polygons in Fig. 3. In Tab. 4 
we report for these two cases the parameters obtained assuming the median values (i.e. 1.0 
and 3.0) of the ranges where their redshift was allowed to vary. 

We can see from Fig. 3 that with the above reasonable assumptions on the uncertainty 
of their spectral parameters or distance also these 8 GRBs are still consistent with the found 
correlation between E pea k and E 7 . 
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5. 



Discussion 



The analysis of a sample of GRBs with measured redshift and peak energy revealed a 
tight correlation between their peak spectral energy and collimation corrected energy: 



The maximum (logarithmic) scatter from this correlation in the sample of 24 GRBs with 
measured jet opening angle is ~0.25 dex while its distribution (Fig. 2(c), insert) has a 
standard deviation of less than 0.1 dex. 

The distribution of the collimation-corrected energy (Fig. 2(b)) for our GRB sample 
presents a wider spread compared to the same distribution in F01 and B03. This is due 
to the different energy bands in which the E 1 has been computed and to the fact that our 
sample includes several recent HETE-2 GRBs not present in the B03 sample (nor obviously 
in the F01 one). Nevertheless the central value of our distribution ~ 6xl0 50 erg is consistent 
with what reported in F01. The tight correlation of Eq. 4 can be found in both the F01 and 
the B03 samples, even though the peak energy range is smaller and therefore the statistical 
significance lower. 

We derived the jet opening angle (Eq. 1) within the framework of the standard afterglow 
theory (Sari et al. 1999). The small scattering around the correlation between E peak and 
E 1 can be seen as an a-posteriori check that the assumptions we made are correct. This 
indicates that the circum-burst density, its radial profile and the efficiency coefficient r^ 7 are 
fairly standard. 

We have checked that no burst among the 8 cases reported in Tab. 3 and Tab. 4, 
for which we have only a partial knowledge of the relevant parameters, contradicts the 
correlation (Fig. 3). One X-ray flash with known redshift (020903) could be consistent 
with our correlation assuming a jet opening angle ~ 25°. The existing radio data are 
controversial and this XRF could be an outlier. The other XRF in our sample (030723) lies 
on our correlation if its redshift is not much smaller than unity. Concerning the two GRBs 
associated with a SN, we have that GRB 980425/SN1998bw (Galama et al. 1998) has a 
very low isotropic energy and no jet break and it is standing alone with respect to any other 
GRB in the E peak -E^,E y>iso plane. On the other hand GRB 030329/SN2003dh (Stanek et 
al. 2003) with a jet opening angle of ~5° (Price et al. 2003a, Tiengo et al. 2003, Berger et 
al. 2004) has a collimation corrected energy and a peak spectral energy which place on the 
Epeak-E^ correlation. 

The sample of 29 GRBs (Tab. 1, 2) also allowed to re-consider the correlation between 
the isotropic equivalent energy E ljiso and the peak energy E peak found by Amati et al. (2002; 




(4) 
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see also Lloyd- Ronning et al. 2000, Lamb et al. 2003; Sakamoto et al. 2004). Including 
all GRBs with measured redshift (Tab. 1) we confirm the existence of such correlation. 
However, we find two differences. On the one hand the spread around the best correlation 
line is larger than previously estimated with smaller samples. On the other hand, the slope 
of the correlation seem to depend on whether XRFs are included or not. We find that 
considering only GRBs, the relation is E pea k oc E®'f®^ om , flatter than previously estimated, 
shown with a dashed line in Fig. f (see also Amati 2004, who finds a slope 0.45 ± 0.06). If 
instead, following Lamb et al. (2003), we include X-ray flashes, the original Amati et al. 
(2002) result fits better the data (dot-dashed line in Fig. 1). The paucity of GRB/XRFs 
with peak energies in the tens of keV range does not allow us to draw any firm conclusion 
(note that the two slopes are consistent at the 2— a level). 

The different slopes of the E peak -E 1 ^ iso (Amati et al. 2002) and the E peak -E y correlations 
imply that they will intersect at some small value of E 1 = E 1 ^ iso (see Fig.l). The precise 
value of the intersection depends mainly on the uncertainties in the E peak -E^^ so relation, as 
illustrated in Fig. 1 (dashed and dot-dashed lines). At any rate, bursts lying in this region 
of the E pea k-E.y ~ E 1 ^ so plane would be characterized, on average, by a nearly isotropic 
emission, small peak energies (possibly in the UV/soft-X-ray range) and by E 7 smaller than 
10 47 erg 3 . 

Recently, it has also been shown (Liang et al. 2004) that the correlation between E peak 
and E 1 ^ iso (found from the time-averaged spectrum) holds when considering the time resolved 
spectral analysis of GRBs. Similarly, due to the much lower scatter (a factor ~ 5.7 lower) 
around the E pea k and E 1 correlation, we expect that considering the time resolved spectral 
peak energy E peak {t) and collimation corrected energy E^it) a possibly even less scattered 
correlation should hold within single GRBs. This, in turn, might help in exploring the (still 
obscure) origin of the prompt emission. In fact, although we do not have an interpretation of 
the found E peak -E 1 correlation, we believe it indicates that there should be some connection 
between the energy emitted by the burst and the emission process which determines its 
spectral properties. 

The existence of the above correlation also allows to predict the jet break time, once 
the spectrum and the redshift of a burst are known. In the forthcoming Swift era this might 
contribute in optimizing the GRB follow up with planned observations at the expected jet 
break time. For this reason we think that it is important that the information about the 
peak energy be disseminated promptly together with the other fundamental GRB properties 



3 Note that GRB 980425 is an outlier whether or not its 7-ray emission was isotropic. 
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such as the position, peak flux etc. Now this routinely occurs for HETE II bursts 4 , and it is 
crucial that this will happen also for the future Swift bursts. 

Berger et al. (2004a) recently discussed the possibility that the real standard energy 
in cosmic explosions (including GRBs, XRFs and hypernovae) is the total released kinetic 
energy rather than the fraction that goes into 7-ray photons. This energy includes corrections 
due to mildly relativistic material (estimated from radio emission such as in SN1998bw; 
Kulkarni et al. 1998) and to a longer time-scale activity of the engine (estimated from 
the afterglow lightcurve, such as in GRB 030329). While not in contrast with this idea, 
our findings underline the fact that the 7-ray energetics is still an important parameter for 
understanding the physics of GRBs. The found correlation implies that the peak energy of 
the spectrum, which is a quantity set by local micro-physical processes, is determined by a 
global property of the jet, namely its total 7-ray energy (or the jet opening angle). 

Last but not least, the E pea k-E 1 relation that we presented in this work and the low 
scatter of the sample of bursts with known redshifts around it strongly reminds a similar 
relation found for SN la between their luminosity and the stretching factor of their optical 
light curve (e.g. Phillips 1993, Perlmutter et al. 1998) with less luminous supernovae showing 
a faster post-maximum light curve decay (Riess et al. 1995). The proper modelling of this 
effect (e.g. Hamuy et al. 1996, Perlmutter et al. 1998) allows to better determine the SN 
luminosity and consequently reduce the scatter in the Hubble diagram with respect to the 
case that assumes that all SN la are standard candles. Similarly, our result shows that among 
GRBs the existence of a correlation with small scatter between their peak spectral energy 
and the collimation corrected energy might be used to determine the different luminosities of 
different bursts. The small scatter around the rest frame E pea k-E-y correlation might be used 
to reduce the scatter of the GRB-Hubble diagram (Bloom et al. 2003, Schaefer et al. 2003) 
and constrain the cosmological parameters (Ghirlanda et al. (2004) in preparation). If our 
relation will be confirmed and possibly extended (particularly to low peak energies) with new 
bursts detected by Swift, we might be able to test the cosmological models farther out the 
predicted SNAP SN limit of z—1.1 (e.g Aldering et al. 2002). This is particularly interesting 
since we might be able to explore an epoch in which the universe, matter dominated, was 
decelerating. 



4 http://space.mit.edu/HETE/Bursts/ 
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6. Conclusions 

We considered all the bursts with measured redshift and spectral peak energy up to now 
(June 2004). We systematically derived their source frame bolometric isotropic equivalent 
energy and, when available from the afterglow modeling (24/38), corrected this energy for 
their collimation angle. We found a very tight correlation between the source frame E peak 
(the peak in the vF v spectrum of the prompt emission) and the collimation corrected emitted 
energy: E peak ~ 480 (£ 7 /10 51 erg) - 7 keV. 

• The small dispersion around this correlation is an indication of the robustness of the 
afterglow theory on which the estimate of the jet opening angle is based (Eq. 1, Sari 
et al. 1999) 

• We found that, differently from previous results, GRBs do not cluster around a unique 
value of their collimation corrected energy, but are spread in a relatively large range 
(~2 orders of magnitudes) centered at E 1 ~6xl0 50 erg. We believe that the smaller 
spread found previously around a similar value ~5xl0 50 erg (Frail et al. 2001, Bloom 
et al. 2003, Panaitescu & Kumar 2001) might be the result of a limited sample of 
GRBs within a relatively small E peak range. 

• The underlying physical motivation of the correlation is mysterious, but it is likely that 
it is the result of the radiation process producing the prompt emission which should 
be related to the local energy content of the burst. 

• There are relatively few bursts with 10 < E peak < 100 keV, just the range of the BAT 
instrument on-board Swift, which therefore will be crucial for confirming or discarding 
the correlation. 

• Finally, and more importantly, the correlation (if the small scatter around it will be 
confirmed) makes GRBs exquisite cosmological tools, to measure Q m , Q\ in a redshift 
range not accessible to Supernovae la, that is, z > 1.7, and up to any redshift. Even 
if the optical information is unavailable for z > 2.5, the necessary information about 
the jet break time can be gathered by following the X-ray and/or the IR afterglow. 
Partial obscuration by dust, which can be an issue for using Supernovae as standard 
candles, does not influence GRBs. 

We thank the referee for her/his comments which improved the paper and Annalisa 
Celotti for useful discussions. DL thanks the Osservatorio Astronomico di Brera for the 
kind hospitality during part of the preparation of this work. G. Ghirlanda thanks the ASI 
(/I/R/390/02) for financial support. We thank M. R. Panzera for technical support. 
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b Redshift determined from: H=Host Galaxy spectrum; XP=X-Ray Photometric data; 0=Optical Transient. 
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c Observed peak energy. 

d References given in order for: redshift, spectral informations. 

References. — (1) Djorgovsky ct al. 1998; (2) Djorgovsky et al. 2001; (3) Kulkarni ct al. 1998; (4) Tinncy 
ct al. 1998; (5) Djorgovsky et al. 1998a; (6)Djorgovsky et al. 1998b; (7) Hjorth 1999; (8) Bloom et al. 2001; 
(9) Vreeswijk et al. 1999; (10) Amati et al. 2000; (11) Galama et al. 1999; (12) Vreeswijk et al. 1999a; (13) 
Andersen et al. 2000; (14) Antonelli et al. 2000; (15) Price et al. 2002; (16) Stanek et al. 2001; (17) Djorgovski 
ct al. 2001a; (18) Infante et al. 2001; (19) Amati 2004; (20) Hjorth et al. 2003; (21) Masetti et al. 2002; (22) 
Price et al. 2002a; (23) Soderberg et al. 2004; (24) Vreeswijk et al. 2003; (25) Greiner et al. 2002; (26) Rol et 
al. 2003; (27) Greiner ct al. 2003; (28) Weidinger et al. 2003; (29) Fynbo et al. 2004; (30) Berger et al. 2004b; 
(31) Amati ct al. 2002; (32)Jimenez et al. 2001; (33) Price et al. 2002; (34) Barraud et al. 2003; (35) Atteia et 
al. 2003; (36) Piro ct al. 2004; (37) Price et al. 2003; (38) Crew et al. 2003; (39) Vanderspek et al. 2004; (40) 
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Table 2. Sample of GRBs 
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a References given in order for: jet break time, external medium density 

b We have assumed that the error is 10% instead of the 1.3% error quoted in Jakobsson et al. 
2003. 

References. — (1) Bloom et al. 2003 (references therein); (2) Bjornsson et al. 2001; (3) Price 
et al. 2002; (4) Jakobsson et al. 2003; (5) Torii et al. 2002, Gorosabel et al. 2002, Bloom et al. 
2002; (6) Delia Valle et al. 2003; (7) Pandey et al. 2004; (8) Andersen et al. 2003; Perterson & 
Price 2003, Burenin et al. 2003; (9)Berger et al. 2004; (10) Tiengo et al. 2003; (11) Bersier et al. 
2004.; (12) Klose et al. 2004; (13) Jakobsson et al. 2004; (14) Dulligan et al. 2003 ; (15) Israel et 
al. 1999 
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Fig. 1. — Rest frame peak energy -E pea fc = E°l s ak (l + z) versus bolometric energy for the 
sample of GRBs with measured redshift reported in Tab. 1. Filled circles: isotropic energy 
corrected for the collimation angle by the factor (1— cos#), for the events for which a jet break 
in the light curve was observed (see Tab. 2). Grey symbols corresponds to lower/upper limits. 
The Solid line represents the best fit to the correlation, i.e. E pea k ~ 480 (£^ 7 /10 51 erg) ' 7 keV. 
Open circles: isotropic equivalent energy E lyiso for the GRBs reported in Tab. 2. The Dashed 
line is the best fit to these points and the dash-dotted line is the correlation reported by 
Amati et al. (2002). 
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Fig. 2. — Top: distribution of the rest frame isotropic equivalent energy E 1}iso for the sample 
of GRBs presented in Tab. 1. Middle: Distribution of the collimation-angle corrected energy 
E y . Bottom: distribution of the distance of each GRB from the correlation between E peak 
and E~f. Insert: distribution of the distance from the correlation with upper and lower limits 
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Table 3. Uncertain cases 



GRB a 


z 


a 


P 


Epeak 

(keV) d 


T90 

(s) 


Fluence 

(erg/cm 2 ) 


Range 
(keV) 


ref c 


970508 S/B 


0.835 


-1.71 (0.1) 


-2.2 (0.25) 


79 (23) 


20 


1.8e-6 (0.3) 


40-700 


(1),(7) 






-1.19 


-1.83 


> 1800 




5.5e-6 [0.5] 


20-2000 


-,(8) 


980326 S/B 




-1.23 (0.21) 


-2.48 (0.31) 


33.8 (17.1) 


9 


0.75e-6 (0.15) 


40-700 


...,(7) 


980329 S/B 




-0.64 (0.14) 


-2.2 (0.8) 


233.7 (37.5) 


25 


6.5e-5 (5) 


40-700 


...,(7) 


991208 I 


0.7 








60 


l.e-4 


25-l.e4 


(2),(9) 


000210 S 


0.846 








20 


6.1e-5 


2-700 


(3),(10) 


000301C I 


2.033 








50 


4.1e-6 


25-l.e4 


(4),(H) 


000418 I 


1.118 








30 


2.e-5 


15-l.e4 


(5),(12) 


000926 I 


2.036 








25 


2.2e-5 


25-l.e4 


(6),(13) 



a Instrument(s) that were triggered by the GRB: S=BeppoSAX; B=BATSE; (B)=also BATSE, 
R=RXTE, A=ASM, I=IPN (Ulysses, NEAR, Konus). 

c References given in order for: redshift, spectral informations 

d Observed Peak energy 



References. — (1) Metzger 1997; (2) Dodonov et al. 1999; (3) Piro et al. 2002; 4) Castro 2000; 
(5) Bloom et al. 2000; (6) Castro 2000a; (7) Amati et al. 2002; (8) Jimenez et al. 2001; (9) Piran 
et al. 2000; (10) Stornelli et al. 2000; (11) Garnavich et al. 2000; (12) Hurley et al. 2000; (13) 
Hurley et al. 2000a. 
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Table 4. Uncertain cases 



GRB a 


tjet 


n 


ref 6 


e 


E ■ 




Epeak(l + z) 




(days) 


(cm -3 ) 






(deg) 


(erg) 


(erg) 


(keV) 


970508 


25 (5) 


3.0 [1-10] 


(1) 




24.0 [3.3] 


7.1e51 (0.15) 


6.1e50 [1.6] 


145 [43] 




25 (5) 


3.0 [1-10] 


(1) 


... 


22.1 [3.05] 


1.4e52 [0.17] 


1.02e51 [0.3] 


>1503 


980326 


<0.4 


3.0 [1-10] 


(1) 




<5.08 


5.6e51 (1.1) 


<2.2e49 


71 [36] 


980329 


<1 


20 [10] 


(1) 


(1) 


<3.33 


2.1e54 (0.2) 


<3.5e51 


935 (150) 


991208 


<2.1 


3.0 [1-10] 


(1) 


... 


<8.64 


l.le53 


<1.2e51 




000210 


>1.7 


3.0 [1-10] 


(1) 




>5.76 


2e53 


>le51 




000301C 


7.3(0.5) 


27(5) 


(1) 


(1) 


13.14 


4.37e52 


1.14e51 




000418 


25 (5) 


27 


(1) 


(1) 


22.3 


7.51e52 


5.6e51 




000926 


1.8(0.1) 


27 (3) 


(1) 


(1) 


6.19 


2.7e53 


1.6e51 





a References given in order for: jet break time and external medium density. 



References. — (1) Bloom et al. 2003. 



-25 - 




Fig. 3. — Rest frame peak energy E peak = E°l s ak (l + z) versus bolometric energy for the 
sample of GRBs reported in Tab. 3,4 and excluded from the larger sample due to (i) uncertain 
redshift (980329, 980326), (ii) uncertain peak energy (970508) or (hi) no published spectral 
information (991208, 000212, 000131C, 000418, 000926). 



